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Abstract 

The de Haas - van Alphen effect in quasi-two-dimensional metals is studied at arbitrary pa- 
rameters. The oscillations of the chemical potential may substantially change the temperature 
dependence of harmonic amplitudes that is usually used to determine the effective electron 
mass. Hence, the processing of the experimental data using the standard Lifshitz-Kosevich 
formula (that assumes the chemical potential to be constant) may lead to substantial errors 
even in the limit of strong harmonic damping. This fact may explain the difference between the 
effective electron masses, determined from the de Haas - van Alphen effect and the cyclotron 
resonance measurements. The oscillations of the chemical potential and the deviations from 
the Lifshitz-Kosevich formula depend on the reservoir density of states, that exists in organic 
metals due to open sheets of Fermi surface. This dependence can be used to determine the 
density of electron states on open sheets of Fermi surface. We present the analytical results of 
the calculations of harmonic amplitudes in some limiting cases that show the importance of the 
oscillations of the chemical potential. The algorithm of the simple numerical calculation of the 
harmonic amplitudes at arbitrary reservoir density of states, arbitrary warping, spin-splitting, 
temperature and Dingle temperature is also described. 

The quantum magnetization oscillations (or the de Haas - van Alphen (dHvA) effect) was 
discovered long agojlj and has been used a lot as a powerful tool for studying the Fermi surfaces 
and single electron properties in metals]^]. In three dimensional (3D) metal the good quantitative 
description of this effect is given by the Lifshitz-Kosevich(L-K) formula || . In two- or quasi- two- 
dimensional compounds the deviations from the L-K formula are possible for three reasons: the 
harmonic damping in two-dimensional (2D) case is different, the impurity scattering may not be 
described by the usual Dingle law and the chemical potential becomes also an oscillating function 
of the magnetic field. The first problem is important only when the harmonic damping is weak and 
can be easily solved by using the 2D harmonic expansion Jl. The second problem concerns with 



an accurate calculation of the density of states (DoS) with electron-electron interactions and the 
impurity scattering. The electron-electron interactions are not very important when many Landau 
levels are occupied (we consider the case when the Fermi energy tp is much greater than the 
Landau level separation and temperature). The impurity scattering in 3D case adds an imaginary 
part *r(i3) to the electron spectrum that means that an electron may leave its quantum state 
with probability w — T(B)/-kT% per second. If one assumes this width T(B) of energy levels to be 
independent of the magnetic field B he gets the Dingle law of harmonic damping^] 

Ai — exp (-2irlT/huj c ) 

where Ai is the amplitude of the harmonic number I and uj c = eB/m*c is the cyclotron frequency. 
This Dingle law has been proved by many experiments on 3D metals. In 2D case this law may 
be incorrect and the problem of the DoS distribution in 2D metals has not been solved yet, 
although many theoretical works have been devoted to this subject (for example, @,f(|, (?]])• The 
problem is complicated because even the exact calculation of the point-like impurity scattering is 
not enough since the long-range impurities (and, probably, the electron-electron interactions) are 
also important in 2D case||. The procedure of extracting the DoS distribution from the dHvA 
measurements has been proposed recently 0. In the present paper we focus on the third question, 
so, we assume the Dingle law to be valid and in this approximation we consider the influence of the 
oscillations of the chemical potential on the harmonic amplitudes of the dHvA oscillations. Since 
we consider the quasi-2D case, the Dingle law is not a bad approximation. We shall show that the 
oscillations of the chemical potential change substantially the temperature and Dingle temperature 
dependence of the harmonic amplitudes even in the limit of strong harmonic damping. Hence, the 
estimate of the effective electron mass based on the L-K formula may lead to the errors up to 30%. 
This can be an explanation of the difference of the effective electron masses obtained from the 
dHvA effect and cyclotron resonance measurements (for example, in JlO[ ] and p|). This problem 
was examined numerically by N. Harrison et al.|12| at zero warping W of Fermi surface(FS). In 
this paper we derive the explicit formulas that describe the quantum magnetization oscillations at 
arbitrary parameters. The analytical study of this result is made in some limiting cases. It shows 
the importance of the oscillations of chemical potential on the harmonic amplitudes. 
The energy spectrum of quasi-two-dimensional electron gas is 

1 W 

E n ,k z ,a = hue (n + -) + — cos(fc 2 rf) + afi e B (1) 

where W is the warping of quasi-cylindrical Fermi surface. The DoS distribution with impurity 
scattering may be written as 

p(E,B)=p (E,B)+p(E,B) 
where the oscillating part of the DoS at E ^> Uuj c isfl3f 



~p(E,B) = P- yVlVcos ( 2ttI-^-] Jo ( nl^-) cos ( 2nl^] exp 

hio c f-' v huj c j \ huj c j \ hw c ) 
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In this formula g = i?/$o is the Landau lcvcl(LL) degeneracy, the factor cos (2irlp, e B /hui c ) is due 
to spin-splitting and the factor Jq (ttIW / ?ko c ) comes from the finite warping W of quasi-cylindrical 
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Fermi surface. Jo (a;) is the Bessel function of zeroth order. The last factor in (g) is the usual 
Dingle factor. 

The non-oscillating part of the DoS 

Po (E,B) = ^-(l + n R (E)) 

where n R (E) is the ratio of the reservoir density of states to the average DoS on quasi-2D part of 
FS. The reservoir density of states exists in quasi-2D organic metals due to open sheets of Fermi 
surface. These quasi-one-dimensional states do not contribute to magnetization oscillations since 
they form a continuous spectrum. 

If the DoS is known one can calculate the thermodynamic potential 



/>oo 

tl(li,B,T) = —T / p(E,B)\n 
Jo 



1 + exp 



dE = a (fx,B,T) + n(fi,B,T) (3) 



T 

where p(B) is the chemical potential and the oscillating part of the thermodynamic potential is|l3[] 

- „ (-1)1 « \l W , ,„ ,u e H. ( 2vr/r\ 

n - 2 5 TV^-^cos (2nl^- — — -Jo (wl— cos (2irl%—) exp — 

f— ^ I TvjJc smn(Ai) nu> c nuj c \ fiui c ) 

where A = 2tt 2 T /Tiu c . The total particle number is usually constant: 

N= f d n(»B,m r p(e,b) E=const 

This is an equation on the chemical potential as a function of magnetic field. Separating the 
oscillating part of the DoS and substituting 

Jo l + exp(^^) 



(ef is the Fermi energy at zero magnetic field) we get 



l>n[ E.B),1E : / P(E ^} N dE (4) 



l + exp(^^) i + exp ^E^y Jo i + exp 



E-jj 
T 



Now we use the fact that the reservoir DoS n R (E) does not change appreciably on the scale of T 
or Ii^ — sfI < hu} c /2 (this is true if many LLs are occupied because tir{E) changes substantially on 
the scale of Fermi energy). Then n R (E) w nji(ep) — const = n R . The left hand side of equation 
([|) now simplifies and after substitution of (||) we get the equation on the oscillating part Ji{B) of 
the chemical potential 

fi{B) = M (5) - e F = ^ c x (5) 
7r(l + n R (eF)) 
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x ±^ (_iy+i ^ / 27r z (£F + A (i?)) \ \i cog / 27r/ ^£\ exp f 27r*r\ Jo / 

^— j' l \ fiLOc J sinh(AZ) \ fiuj c ) \ Uu) c J \ fou> c J 

This nonlinear equation can not be solved analytically without any approximations but it de- 
termines the oscillations of the chemical potential arbitrary parameters (only ef S> T,Tilu c is 
assumed) . 

The magnetization oscillations at constant electron density N — const 



M 



d{n + N[i) _ on 



N— const o j-y \fi.N— const 



dB dB 



q \N,B— const "T-tV I \N— const QB Imj-^— cons * 

The oscillating part of the magnetization 

M ^ = ~ m U=~n.,= £e F E / s . nhA? cos(2^)exp (- — j x (6) 

nuj c nuj c 2/i nw c nxj c 

where is given by equation (||) and contains the dependence of magnetization on the reservoir 
DoS. The formulas (||) and (^) describe the magnetization oscillations at arbitrary parameters. The 
only approximation, used in these formulas is the Dingle law of harmonic damping. In quasi-2D 
organic metals with warping W > To the Dingle law is believed to be a quite good approximation. 

The formulas @) and (||) are the good starting point for the numerical calculations. From these 
formulas we see that in the limit W//J, <C 1 the oscillating parts of magnetization and chemical 
potential are connected by the simple relation 

At zero warping this was obtained in || . 

The nonlinear equation (g) for Ji(B) can be solved analytically only in some simple approxi- 
mations. We shall do this to illustrate the influence of the oscillations of the chemical potential on 
the temperature and Dingle temperature dependence of the harmonic amplitudes. So, we consider 
zero warping, zero spin-splitting and zero temperature. Then the sum in the right-hand side of 
equation (||) can be calculated and we get 

x 1 ( sinfy + x) \ . . 

■ arctan , x », (7) 



2 (1 + ur) \cos(y + x) + e 

where x = 2irji(B)/hLo c , y = 2-KEp /Hu> c and b = 2ttT/?uv c . 

At very large electron reservoir tir = oo , x = and we have the limit of fixed chemical 
potential. In this case the magnetization is given by E3] 

M(B) = ^ arctan ( b ^ v \ - ) (8) 
V ' ttB \e b + cos(y) J y ' 
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The temperature dependence of the harmonic amplitudes is given by the L-K formula: 

2^Tl/f^ c 
lK ' sinh(27r 2 TZ/?iw c ) 

It is possible to solve analytically the equation (Q) also at ur = and ur = 1. At zero electron 
reservoir ur = the solution of this equation is 

x It-(B) ( sin(y) 

- 7r -7- — - = arctan ' 



2 hu> c \e b — cos(y) / 

It gives the oscillations of the chemical potential. The magnetization at zero electron reservoir 

2ge F ( sin(y) 

arctan —r ; . 

nB \e° — cos(y) 



M(B) = arctan ( ^ m ^ - ) (10) 



It coincides with (||) after the phase shift y — * y + tt and the sign change M — > — M . This means 
that the harmonic damping law 

Ai ~ (I//) • exp (-/ • 6) (11) 

does not change, and only the sign of all even harmonics turns to inverse. This symmetry between 
the cases of fixed chemical potential fi = const and constant particle density N = const is a feature 
of special exponential law of harmonic damping. Any finite temperature and the density of electron 
reservoir breaks this symmetry. 

Let us consider now the intermediate case ur = 1. The equation (fjj) becomes 

sin 2 sm(y + x) 



cosx cos(y + x) + e b 



It gives 



arcsm 



in (e siny 



and the magnetization becomes 

M(y) = ^—^ arcsin (e^siny) (13) 

To say how the harmonic damping has changed we have to calculate the amplitudes of several first 
harmonics of this expression. The amplitude of the first harmonic is 



1 f n 

AAb) = — I arcsin (e _b sin y) sin ydy 
after integration by parts we get 

4 r /2 cos 2 ye~ b dy 



4 r 
Ai (b ) = - \ 

I" JO 



VT— e" 2b sin 2 y 
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This is the superposition of two elliptic integrals: 



Ax (6) = - [e b E(e~ b ) - 2 sinh b K( e - b )] (14) 



At b ^> 1 the deviations of A\(b) from the L-K formula are small: 

A 1 (b)=e- b + e- 3b /8 + .. 



In the opposite limit b <C 1 we get 



Al (b) = ~\l-b(\n4=-\)+0{tf)} ! 1 ■") ! 



7T I V V26 2 

This is substantially different from the L-K dependence Ai(b) — exp(— b) m 1 — b. For example, 
the value A 1 (Q) is 4/ir times larger than the L-K prediction. 



The amplitudes of the next harmonics reveal the stronger deviation from the L-K formula (11). 
All even harmonics disappear since the expression ( fl3| ) has the symmetries M(ir — y) = M{y) and 
M(-y) = M(y). 

The amplitude of third harmonic can also be calculated. At b ^> 1, e~ b <C 1, 

A 3 (b) = -fr 3 7l2 + 0(e- 5& ) 

This result is in contrast to the cases ur = or tir = oo where we had A3 (b) — e _3b /3 . It is 
not surprising since in the symmetric case ur — 1 the oscillations should be much smoother and 
more sinusoidal. Hence one should have an increase of first harmonic and a decrease of higher 
harmonics. At b = 



4 W2 

A 3 (0) = - / 



cos 3;/ cos ydy = 4_ 

cos y 9n 



which is ~ 2.35 times less than the Lifshitz-Kosevich prediction ^3(0) = 1/3 and has the opposite 
sign. So, in the case ur = 1 the first harmonic is increased while the others are strongly decreased 
in amplitude compared to the cases of zero and infinite electron reservoir. The deviation from 
the L-K formula reduces as the warping of FS increases. The above analysis shows also that the 
harmonic ratios at low temperature and Dingle temperature can give a quantitative estimate of 
the electron reservoir density which is much more precise than just a note about the slope of the 
magnetization curve. 

To include correct temperature dependence, warping and spin-splitting and consider arbitrary 
reservoir density one can perform the numerical calculations, based on the solution of the equation 
(||) for the chemical potential and substitution of this solution into the formula (||) for the magne- 
tization. This can be done at arbitrary parameters, available on the experiment. The temperature 
dependence of the first three harmonic amplitudes is given in fig. 1 for the following set of param- 
eters, close to the real experiments on a-(BEDT — TTF)2KH g{SCN)i\ the reservoir density 
ur = 1, the dHvA frequency F = 700 Tcsla, the effective mass m* = 2?7io, the Dingle temperature 
Td = 0.2K and the warping W = IK. We see the substantial deviation from the Lifshitz-Kosevich 
dependence. The obtained amplitude of the first harmonic at T — * is about 1.1 times larger 
than the Lifshitz-Kosevich prediction. If we put also Tjj — ► and W — > their ratio becomes 
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Figure 1: Temperature dependence of harmonic amplitudes. The solid lines are the numerical 
results (at tir = 1, m* = 2mo, Td — 0.2K and W = IK; see text) and the dashed lines represent 
the Lifshitz-Kosevich prediction at the same parameters. Their strong deviations are clearly seen, 
especially for higher harmonics. 

4/7T = 1.27 in agreement with the analytical result (|l5|). The second harmonic amplitude is close 
to zero at T = 0. The amplitude of the third harmonic changes the sign at T w 0.8K and deviates 
very strongly from the L-K formula. It is damped much stronger than the L-K predictions. At 
T = and W = it also coincides with the prediction (|l6|). 

To conclude, it was shown both analytically and numerically that the oscillations of the chemical 
potential are essential for the temperature dependence of harmonic amplitudes of dHvA oscillations 
in quasi-two-dimensional compounds. The accurate determination of the effective electron mass 
from the dHvA effect should take this effect into account. This can be done by the simple numerical 
calculation based on the formulas (||) and (0). The oscillations of the chemical potential depend on 
the reservoir density of states according to the formula (||). This fact may be used for the estimate 
of the reservoir density of states in organic metals. 

The author thanks A.M. Dyugaev and M.V. Kartsovnik for useful discussions. The work was 
supported by RFBR grant N 00-02-17729a. 
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